Superfluidity of Dense 4 He in Vycor 
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We calculate properties of a model of 4 He in Vycor using the Path Integral Monte Carlo method. 
We find that 4 He forms a distinct layered structure with a highly localized first layer, a disordered 
second layer with some atoms delocalized and able to give rise to the observed superfluid response, 
and higher layers nearly perfect crystals. The addition of a single 3 He atom was enough to bring 
down the total superfluidity by blocking the exchange in the second layer. Our results are consistent 
with the persistent liquid layer model to explain the observations. Such a model may be relevant to 
the experiments on bulk solid 4 He, if there is a fine network of grain boundaries in those systems. 
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There is a long history of experiments of helium ab- 
sorbed in porous media carried out to probe the response 
of superfluidity to disorder |lj . Recently, Kim and Chan 
reported the observation of a supersolid Helium phase, 
in porous Vycor Q and in bulk solid 4 He Q. A su- 
persolid is a proposed phase of a quantum system in 
which long-range crystalline order and superfluidity co- 
exist. We focus here on the measurements of the Helium- 
Vycor system. 

The occurrence of supersolid behavior in Vycor, a dis- 
ordered porous glass, could be understood from the prop- 
erties of Vycor. One explanation is that the complex 
Vycor geometry stabilizes mobile defects which then un- 
dergoe BEC at low temperatures « 0.2K. However, since 
the phenomenon is observed in bulk solid 4 He it is not the 
Vycor pores that are essential. Although, Kim and Chan 
pressurized their cell to 60 Bars, substantially above the 
estimated freezing pressure for helium in Vycor of 40 
bars, even if most of the helium is solid, it is not clear if 
there remains a liquid film (the persistent liquid layer or 
PLL) near the Vycor-helium surface. The film could arise 
from the mismatch of lattice parameters as the density of 
4 He varies from the center of the pores to their surface. 

Even though there have been many experimental stud- 
ies of helium in Vycor, there have been few microscopic 
calculations of the detailed microscopic structure of this 
system. Here we report on a model of the helium- Vycor 
system and calculate its properties with the Path Inte- 
gral Monte Carlo(PIMC) method 0. PIMC can calcu- 
late exact thermodynamic properties of bosonic systems 
such as He at non-zero temperature by sampling the 
thermal density matrix p = e~^ R , with (3 = l/kpT and 
H the Hamiltonian. An explicit expression for the den- 
sity matrix is obtained by expanding into a path and 
approximating the higher temperature density matrices. 
Bose statistics are obtained by symmetrizing the density 
matrix p Bose {R, R' ; (3) = J2p p( R , PR 'i/ 3 ) l Nl PIMC 
proved accurate in studying properties in the normal 
liquid, superfluid and crystal phase [jj. In contrast to 
methods based on trial wavefunctions, in PIMC only the 
Hamiltonian enters, so no biases are made on the struc- 



ture of the many-body system. 
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FIG. 1: The external potential V(z) experienced by the he- 
lium atoms in "cell 2". The Vycor is on the left side (z=- 
12. 5A) and solid helium on the right side (z=12.5A). The 
inset is a 3D representation of the rough Vycor surface: the 
black dots are the positions of the Vycor impurities placed 
randomly at lA away from the Vycor wall. The rugged sur- 
face shows the positions of the helium atoms (located at the 
vertices) in the first layer. 

Vycor's complex geometry is difficult to simulate di- 
rectly. Under the transmission electron microscope, Vy- 
cor shows interconnected pores with diameters between 
60A and 80A and a length of 300A, with a narrow dis- 
tribution of pore sizes. Current simulation techniques 
treating all the atoms with PIMC are unable to simulate 
even a single pore (containing roughly 30000 atoms). It 
is thought that above the freezing pressure (40 bars) a 
4 He crystal will occupy the center of a pore Previ- 
ous PIMC calculations have found that a perfect crystal 
will not have a superfluid response at long wavelength 
7]. Making this assumption, we study only the region 
near the surface of a pore and model it with the slab 
geometry: our simulation cell is periodic in both the x 
and y directions. In the negative z-direction there is a 
wall, representing bulk Vycor on top of which are placed 



Vycor "particles" in random positions 1 A above the flat 
Vycor surface. These particles serve to make the surface 
rough and break translational invariance. In the positive 
z-direction is a wall which models bulk solid Helium at a 
pressure of 62 bars (see Fig.QJ . 

The main input to PIMC is the interaction potential 
between various particles. We assume that the Helium- 
Helium interaction is given by the Aziz 8] potential 
which has been used within PIMC to study numerous 
other properties of helium and gives energies accurate 
to about 1% of the bulk binding energy. The potential 
between helium and the upper wall was derived by inte- 
grating an approximate LJ 6-12 potential (e = 10.22K 
and a = 2.556A) over the volume z > 11. 7 A (cell 1) 
resulting in a LJ (3-9) potential. 

We also assume helium- Vycor surface interaction is a 
LJ (3-9) potential 0: V{z) = f [(f) 9 - 3(f) 3 ]. Since 
Vycor glass, SiC>2, should behave similarly to MgO [loj . 
we pick the well depth to be D = — 86.9K, and the range 
of the attraction, z e = 3.6A In order to pin the helium 
crystal in the xy plane and model the roughness of the 
Vycor, we add Vycor "particles" in random positions 1 
A above the wall. The interaction between the helium 
atoms and the Vycor particles is determined by demand- 
ing that a complete layer of the particles give the same 
LJ (3-9)potential. Experiments [ll| show that the rough- 
ness in Vycor is on the scale of 0.8nm, not very different 
from what we have assumed. Figure ^ shows the wall 
potential and the surface roughness. 

We set the helium density to match the experimental 
conditions by adjusting the number of helium atoms and 
the total area in the xy plane so that the helium den- 
sity in the topmost layer matches that of solid 4 He at a 
target pressure close to 62 bars. We start the simulation 
with 221 atoms placed in 7 layers in an hexagonal closed 
packed solid phase. Each layer contains 30 atoms except 
for the first layer placed at the strongly attractive Vy- 
cor wall potential (Fig.QJ. We have performed extensive 
simulations with two geometries, denoted as Cell 1 (221 
He atoms with a box 17.75 x 18.45 x 23.4A) and Cell 2 
(221 atoms with a box 17.25 x 17.93 x 25A). Cell 1 is 
roughly stress free, while Cell 2 provides us with a way 
to look at the model under anisotropic stress. 

The density in the z-direction (see Fig.[2J shows a dis- 
tinct layered structure of 4 He. The density increases as 
we approach the Vycor wall because of the stronger at- 
traction of the potential well. The Vycor particles distort 
the shape of the density peak in the first layer because 
there are binding sites at different values of z. We find a 
perfect crystal is stable in the upper portion of the cell 
with a lattice constant of 3.55 A in cell 1 and 3.45 A in 
cell 2, at the target pressure of 62 bars. We estimate the 
pressure by calculating the density per layer from fig. [21 
and comparing it with the experimental equation of state 

We determine the spatial ordering within a layer with 
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FIG. 2: The density and superfluid density as a function of 
z. The left scale (solid line) corresponds to pure 4 He and the 
right to the mixture: one 4 He atom replaced by a 3 He atom. 
Top: The density in the z-direction for cell 1 at 0.2K. The split 
peak in the first (leftmost) layer is due to the rough Vycor 
surface. The dotted lines show the 3 He density at 0.2K and 
0.7K. Bottom: Local superfluid density of 4 He determined by 
recording which layers the winding paths visit; Eq. (3). 

the structure factor in the x-y direction 

Sn(k) = l( Pn (k) Pn (-k)) (1) 

where j0 n (k) = Y^H=x®{ z i e n)exp(ik • r^) is the 
Fourier transform of the density within layer n and 
k = (k x , k y ,0). We can see the signature of a solid from 
the peak of S(k) around fco=2.04A _1 as shown in fig. 01 
for cell 1; the peak clearly shows an hexagonal structure 
for layers three and above. The density profiles within a 
layer at the top of the figure confirm this interpretation. 
However, not all of the layers are solid. In fact, layer 1 is 
solid-like with the helium atoms well-localized but with 
their mean positions determined by the underlying dis- 
order. Layer 2 is more disordered, and the atoms are out 
of register with the first layer. Because the second layer 
density is lower, the atoms are much less localized and, 
as we shall see, are able to become superfluid. Layers 
three and above are quantum solids, mostly free of de- 
fects. This is relatively independent of pressure, because 
of the strong Vycor potential shown in Fig. ( 1). We note 
that recent neutron scattering experiments |14j support 
the co-existence of solid-like and liquid- like layers at these 
pressures. 

The superfluid fraction is computed in PIMC from 
the mean squared winding number |l3j |. p s /p = 
m(W 2 )/(2ph 2 N) where N is the number of Helium 
atoms and W = 2ik( ri > k — r i,k+i); the sum is over 
particles i, and time slices k. The superfluid fraction in- 
creases as we lower the temperature as shown in fig. 0] 
and approaches values of about 4% below 0.3 K. 

To find the spatial distribution of the superfluid den- 
sity, we divide the winding number estimator into local 
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FIG. 3: Top: Contour plots of the helium density in the x-y 
plane in the first three layers at 0.2K for cell 1. In the first 
layer above the Vycor, the atoms are pinned by the strong 
Vycor interaction, in the second layer they are relatively delo- 
calized, whereas the higher levels show a density distribution 
characteristic of bulk solid 4 He. Bottom: The layer-specific 
structure factor. The first layer has an amorphous structure, 
the second is still distorted but with a much smaller peak, 
while the third and higher layers have a single large peak at 
k o =2.04A _1 characteristic of a 2D quantum solid. 



contributions that sum to the total superfluid density 
(lif . The superfluid density as a function of the distance 
above the Vycor wall x( z ) i s: 
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In fig. we show X {z) for T=0.2K in cell 1. One 
can see the layered structure of the density. Layer 2 has 
the largest superfluid component. Layer one contributes 
because some atoms sit close to layer 2. Layer 3 is also 
active. The superfluid response goes to zero above layer 
3. However, this decay is slower in cell 2 where the two 
additional layers (4 and 5) still contribute to p s /p- 

To compare with experiment we must make two correc- 
tions; first that our model has a larger fraction of atoms 
near the Vycor surface than experiment, and second that 
our cell has no tortuosity: the experimental path length 
for a superflow is greater than the straight line distance. 
We correct for the difference in Vycor surface area by 
assuming that the superfluidity is confined to a finite 
distance above the Vycor so that the effective number 
of atoms that do not respond to the moving boundaries 
is: N s = aA where A is the Vycor surface area (actu- 
ally the surface area of the mobile layer) . We find that 
the maximum number of superfluid atoms per unit area 
is: a w 2.7nm~ 2 . Then for a Vycor sample with exper- 
imental surface area per unit volume of a = 0.2nro -1 
the relative superfluid response would be aa s //j,pK where 
p = 0.3 is the experimental Vycor sample pore fraction, 
k = 5 is the tortuosity, and p = 31.66nm -3 is the solid 
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FIG. 4: The superfluid fraction ps/p vs. temperature. The 
upper curve is for pure 4 He. The lower curve is for the system 
with a 4 He atom replaced by a 3 He atom. 



number density. Using the experimental parameters of 
Vycor, we would predict a measured superfluid fraction 
of 0.011. In fact, Kim and Chan measure a value two 
times smaller than this. One important effect missing 
in our calculation is that the supercurrents must tunnel 
through various weak links, an aspect not present in our 
calculations because of the limited extent of the cell in 
the x-y directions. In addition, we have not taken the 
thermodynamic limit, though it is not obvious how to do 
this without a better model of the Vycor. We also note 
that the model is lacking important aspects of disorder, 
for example, that induced by the curvature of the pore 
and possible grain boundaries between the pores. Given 
these uncertainties, Kim and Chan's measured values are 
well within the range expected from our calculations. 

Recent experiments [l6| do not find any indication that 
the liquid layer can flow in response to a pressure gra- 
dient. It is possible that there is an "insulating" {e.g. 
defect-free, see below)layer of solid 4 He which forms at 
the surface of the Vycor samples and blocks the super- 
flow. 

Kim and Chan also studied the effect of 3 He impurities 
on the measurement of the superfluid density. Accord- 
ingly, we replaced a single 4 He with a 3 He atom. (Note 
that we do not need to consider fermi statistics for a 
single fermion). This corresponds to roughly 0.3% con- 
centration of 3 He (assuming our cell size is ks 2/3 of the 
pore) . Experimentally, this concentration was enough to 
destroy any "supersolid" response. With local moves the 
system was slow to reach equilibrium: we used a move 
that swapped the identity of the impurity with a ran- 
dom 4 He atom. PIMC with the new procedure quickly 
reached equilibrium. 

The impurity, being lighter than 4 He, has a larger ki- 
netic energy and hence a larger atomic volume. This 
extra space is available in the second liquid-like layer. 
Also, the 3 He atom does not exchange with the 4 He and 
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thus reduces the superfluid density in its neighborhood. 
The density profiles show that the 3 He impurity prefer- 
entially goes to the same site where the superfluidity is 
maximized (Fig. [2Jl . At higher temperatures, it tunnels 
to other layers with an excitation energy of 0.8K. As the 
temperature is lowered, two competing effects take place. 
The 4 He atoms nearest to the Vycor surface exchange re- 
sulting in superfluidity. However, the 3 He atom migrates 
towards the most superfluid layers, hence reducing the 
superfluid response and a shift towards a lower transi- 
tion temperature, see fig. 21 It is plausible that in a 
larger cell, an even smaller concentration of 3 He could 
pinch off the winding exchanges (or supercurrent) by go- 
ing to choke positions not present in our small cell, thus 
giving better agreement with the experimental findings 
of a critical 3 He concentration of 0.1%. 

Kim and Chan mention two pieces of evidence to 
argue against the liquid film interpretation. First, the ob- 
served temperature dependance of the superfluid density 
is unlike that of films. However, the films under pressure 
are totally enclosed within a solid and they are not like 
films forming an interface between a solid and a vacuum; 
as we have seen, there are low energy excitations giving 
rise to pronounced temperature effects, not present in 
the later case. Also, the connectivity whether primarily 
2D or 3D, of these solid-enclosed films could be different; 
it is likely that they are gossamery, as opposed to the 
robust films resulting when the pores are only partially 
filled. The second effect mentioned by Kim and Chan is 
the poisoning by small amounts of 3 He. This does not 
happen in free films because the 3 He atom is above the 
plane of the 4 He film in a delocalized state and thus is 
not effective in preventing 4 He exchanges. 

The liquid-layer scenario can be used to give insight 
into the other porous media and bulk helium experi- 
ments. Let us assume that the observed superfluidity 
is due to surface superfluidity as proposed by Burovski 
at al. 0| and that the helium microcrystals are roughly 
spherical with a diameter of _R, giving a surface to volume 
ratio of a = 3/R. Arguing as before, to obtain the ob- 
served 2% superfluid response of bulk 4 He, we must have 
that R = 3cr/(Kp s ) where a is the superfluid density per 
unit area at the crystal interfaces. It is reasonable to ex- 
pect that a is considerably larger than at the interface 
between helium and Vycor, since the pressure exerted by 
the silica increases the density and hence decreases the 
mobility of the helium atoms. Let us assume that there 
are 2 complete superfluid layers at the interface giving 
a fa 20nm~ 2 . Given that the grain boundaries are larger 
and less fractal, we might expect that k rj 2. Using the 
experimental p s for bulk He, we find the average grain 
size, R — 50nm. It is quite likely that 3 He will stabilize 
the grain boundaries, thus explaining how such a small 
concentration of 3 He can affect the response. Such spec- 
ulation need to be confirmed by performing experiments 



on much better crystals and measuring the density and 
the sizes of grains. 

Our results show superfluidity is localized in specific 
layers of 4 He above a Vycor surface. We obtain a super- 
fluid response about 2 times what is observed, but the 
difference is likely due to the very small simulation cells 
we used that do not have the full range of the random 
disorder and crystal defects. We also found that 3 He 
impurities gravitate to the same spatial locations as the 
superfluid density, thus poisoning the effect. Based on 
these simulations, the persistent liquid layer interpreta- 
tion of the Kim-Chan experiment seems not to be ruled 
out. Further studies with larger cells and more realistic 
disorder are needed to firm up these conclusions and to 
see if this mechanism, when involving grain boundaries, 
could give rise to the phenomena in bulk solid helium. 
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